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Abstract 

Interactions of noncommutative solitons in a modified U(n) sigma model in 2+1 dimensions 
can be analyzed exactly. Using an extension of the dressing method, we construct explicit time- 
dependent solutions of its noncommutative field equation by iteratively solving linear equations. 
The approach is illustrated by presenting bound states and right-angle scattering configurations 
for two noncommutative solitons. 
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1 Introduction 



Noncommutativity of spatial coordinates offers a way to introduce nonlocality into (quantum) field 
theory without losing control over its structure. Before attempting to quantize noncommutative 
field theories, it is certainly warranted to characterize the moduli space of their classical configura- 
tions. The generalization of solitons and instantons to the noncommutative realm has just begun. 
For recent reviews on the subject and its connection with string and brane theory see [|l], & [|, [I| 
which include the original literature. 

In H], ^| we have shown that open N=2 strings in a l?-field background induce on the world 
volume of n coincident D2-branes a noncommutative generalization of a modified U(n) sigma model. 
The topological nature of N=2 strings and the integrability of their tree-level dynamics render this 
noncommutative sigma model integrable. In || we have described a family of multi-soliton 
solutions to its noncommutative field equation. These multi-soliton configurations, explicitly given 
for finite values of the noncommutativity parameter 9 and for any U(n) gauge group, represent 
q lumps of energy moving with constant velocities and escaping completely unharmed from their 
mutual encounters, i.e. they do not interact. 

The question arises whether the no-scattering feature of these solitons is a consequence of the 
integrability of the theory or merely due to the restrictions we imposed on the form of the solutions. 
This paper affirms the second explanation: We show how to explicitly construct noncommutative 
multi-soliton configurations with nontrivial interaction within the modified sigma model, general- 
izing the class of solutions presented in || Q • 

The no-scattering configurations were obtained with a (noncommutative variant of) a solution- 
generating technique called the 'dressing method' [||, ||, |l0|, [II], |j~2fl , by using for the tp function 
an ansatz with only first-order poles in the spectral parameter (. Here, our goal is to extend this 
approach to the case where ip possesses higher-order poles. We shall see that this property leads 
to genuine soliton-soliton interaction. 

After presenting the model and our conventions concerning noncommutativity, we briefly review 
the dressing approach, specialize it to the cases of interest, and reduce it to a couple of linear 
equations. As examples we construct two kinds of (abelian and nonabelian) solutions: (a) a two- 
soliton bound state exhibiting a time-dependent ring-like shape and (b) two solitons colliding 
head-on and scattering at 90°. To our knowledge these are the first exact finite-^ multi-soliton 
scattering configurations in 2+1 dimensions. 

2 Noncommutative modified sigma model in 2+1 dimensions 

Commutative model. As has been known for some time, nonlinear sigma models in 2+1 di- 
mensions may be Lorentz-invariant or integrable but not both. In this paper we choose the second 
property and investigate the noncommutative extension of a modified U(n) sigma model (so as 
to be integrable) introduced by Ward |13[| . This model describes the dynamics of a [/(re)-valued 
field <I>(i, x,y) whose arguments are the (real) coordinates (t,x,y) of R ,2 with Minkowski met- 
ric diag(-l, +1, +1). 

The commutative field equation for this modified sigma model reads 

-d t ($- 1 ^$) + d x {$- l d x $) + d y (s- 1 ^) + d y {<s>- x d t $) - d t = o , (i) 

the last two terms of which explicitly break the Lorentz group of SO(l, 2) to the GL(1, R) generated 
by the boost in y direction. Nevertheless, the model possesses a conserved and rotationally invariant 
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energy functional, 

E = JdxdyS = ^jdxdy tr (d t &d t $ + d x &d x $ + d y &d y <f>} , (2) 

where 'tr' implies the trace over the U(n) group space, and <J>t = 

The field equation (|l]) can be obtained from the self-dual Yang-Mills equations in R 2 ' 2 by 
dimensional reduction and a (non-covariant) choice of gauge. In this interpretation, $ arises as 
a (Yang- type) prepotential ]l4| for the Bogomolnyi equations of the 2+1 dimensional Yang-Mills- 
Higgs system 0. 

Noncommutative model. The noncommutative extension of this classical field theory is achieved 
by deforming the ordinary product of classical fields (or their components) to the noncommutative 
star product 

(f*g)(t,x,y) = f(t,x,y) exp{-(d x 9 d y - d y 9 d x }g(t,x,y) , (3) 

where 9 is a positive real constant. Note that the time coordinate remains commutative. In the 
convenient coordinate combinations 

u := \{t + y) , v := \{t - y) , d u = d t + d y , d v = d t - d y (4) 

the noncommutative field equation becomes 

dx^-Ud^-dvi^-Udu®) = o . (5) 

The nonlocality of the star product renders explicit computations cumbersome. We therefore 
pass to the operator formalism, which trades the star product for operator-valued spatial coordi- 
nates (x,y) or their complex combinations (z,z) = (x + iy,x — iy), subject to 

M] = [t,y] = , [x,y] = W => [z,% = 29 . (6) 

The latter equation suggests the introduction of creation and annihilation operators, 

z and a) = ~z with [a, a)} = 1 , (7) 



29 V29 

which act on a harmonic-oscillator Fock space TC with an orthonormal basis {\n), n = 0, 1,2, . . . } 
such that 

a) a\n) =: N \n) = n\n) , a\n) = \fn\n— 1) , a'|n) = y/n+1 |n+l) . (8) 

Moyal-Weyl map. Any function f(t,z,z) can be related to an operator-valued function f(t) = 
F(t,a,aJ) acting in Tt, with the help of the Moyal-Weyl map 

f(t,z,z) — ► F{t,a,a ] ) = Weyl-ordered f{t,VWa, \f29 a ] ) . (9) 

The inverse transformation recovers the c-number function, 

F(t,a,<J) — > f(t,z,z) = F,(t,^=,^=) , (10) 
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where F+ is obtained from F by replacing ordinary with star products. Under the Moyal-Weyl 
map, we have 

dzf — dj = ^= [at,/] and d,f — > 9 y /=-i=[a,/] , (11) 



/*« — » /$ and /dxdy/ = 2vr^Tr/ = 2tt0 £<n|/|n> , (12) 

J n>0 

where 'TV' signifies the trace over the Fock space 7i. For notational simplicity we will from now on 
omit the hats over the operators except when confusion may arise. 

3 Dressing approach and explicit solutions 

The payoff for considering an integrable model is the availability of powerful techniques for con- 
structing solutions to the equation of motion. One of these tools is the so-called 'dressing method', 



which was invented to generate solutions for commutative integrable systems [||, 10, O, |12| and 
is easily extended to the noncommutative setup || [?J . Let us briefly present this method (already 
in the noncommutative context) before applying it to the modified sigma model. 

Linear system. The key observation is that our field equation can be obtained as a compatibility 
condition for a linear system. We consider the two linear equations 

(Cd x -d u )i{j = Aifj and (£d v - d x )i{j = , (13) 

where ip depends on (t, x, y, C) or, equivalently, on (x, u, v, £) and is an nxn matrix whose elements 
act as operators in the Fock space TC. The matrices A and B are of the same type as V but do 
not depend on £. The spectral parameter £ lies in the extended complex plane. The matrix ip is 
subject to the following reality condition |y 



^(t,x,y,C) mt,x,y,()]^ = 1 , (14) 

where 'f ' is hermitean conjugation. The compatibility conditions for the linear system of differential 
equations (|i~3|) read 



d x B-d v A = , (15) 
d x A - d u B - [A, B] = . (16) 

We can solve the second equation by putting 

A = and B = d x ® , (17) 

which transforms the first equation into (the operator version of) our Yang- type equation (||), 

d x i^ 1 a x $) - d v c*- 1 a u $) = o . (is) 

Inserting the parametrization of A and B into the linear system (^) we immediately obtain the 
standard conditions [15] 



^(t,x,y,C^oo) = 1 + OiC 1 ) , (19) 
^(t,x,y,(^0) = <S>-Ht,x,y) + 0(0 . (20) 
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Note that the second equation yields $ directly in terms of tp and thus also A and B via (|17|) or 
directly from 

-i/}(t,x,y,0 ((d x -d u )[tl;(t,x,y,Oy = A(t,x,y) , (21) 
-iP(t,x,y,() {Cd v -d x )m,x,y,C)\ ] = B(t,x,y) . (22) 

Ansatz. Having identified auxiliary linear first-order differential equations pertaining to our 
second-order nonlinear equation, we set out to solve the former. In two previous papers ||. [?J 
we considered solutions ip of the linear system (13) containing only first-order poles in £, 

^ = 1 + tAz > (23) 
k=i ^ 

where the pLf. are complex constants and the nxn matrices Rk(t,x,y) are independent of £. To 
such ip there correspond solutions A and B of equations ( |i~5|) and ([!(]) which are parametrized by 
some n x r matrices T^{t, x,y) || [7j. 

The simplest case occurs when tfj has only one pole at £ = — i, 

-0 = 1 + — =: 1 — P so that $ = 1 - 2 P . (24) 

C+i C+i 

In this case all configurations are static and parametrized by a hermitean projector P = T^J^T^ 
which obeys 

(l-P)8 g P = (l-P)aP = . (25) 

Then T satisfies the equation 

(l-P)oT = , (26) 

which means that aT lies in the kernel of 1— P. Recall that T is an n x r matrix. In the abelian 
case n = 1, and r is the rank of the projector P in the Fock space 7i. In the nonabelian case n > 2, 
and r(< n) can be identified with the rank of the projector in the U(n) group space (but not in 
TL). In terms of P the matrices A and B from the linear system (^) are expressed as 

A = -2id x P and B = 2id y P . (27) 

Dressing approach. The dressing method is a recursive procedure for generating a new solution 
(ip,A,B) from an old one, (ip,A,B) by multiplication, 

m R 

i>{t,x,y,() = x{t,x,y,()ii(t,x,y,() with % = 1 + 7Z ' ( 28 ) 

k=i ^ 

We now choose a static configuration ijj as a seed solution and consider a dressing transformation 
with x being of the same form as tp, 

^-^ = ^ = (1- Ap) (1 _ Ji-P) = i _ ^t(P + P) - J£^ji pp > W 
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where P is some matrix. This leads to a configuration of the form 

$ = ^" 1 (C=0) = (1 - 2P) (1 - 2P) . (30) 

Since both ip and x possess a first-order pole at £ = — i our ansatz for ip contains a second-order 
pole. Substituting ( |29| ) into the reality condition (jlj) directly yields 

Pt = P and P 2 = P , (31) 
qualifying P as a hermitean projector, i.e. 

P =f-J-.T t (32) 
ptp 

with some n x f matrix T. 

Linear equations. Demanding that ^ is again a solution of the linear equations ( |l3| ) with some 
A and P, we obtain 

A{t,x,y) = -ip(t,x,y,() {C,d x -d u ) [ip(t, x, y, C)] 1 " 

= (l-Ap) i( i + JLp) _ (i_ * p)( Caa! _0 u )(i + JLp) , (33) 

C+i C-i C+i C-i 

B(t,x,y) = -4>(t,x,y,() {C,d v - d x ) $(t,x,y,£)]i 

= (^-^P)B(l + ^P) - (l-^P){Qd v -d x ){l + ^P) . (34) 

The poles at C, = ±i on the r.h.s. of these equations have to be removable since A and B are 
independent of £. Putting to zero the corresponding residues, we obtain the essential equations 

(l-P){d g P + (d,P)P\ = and (1-P) {\d t P + {d z P)P} = . (35) 
With the help of the identities 

(l-P)P = and (l-P)f = (36) 
the equations ( |35| ) reduce to 

(1-P) [at + [a, P] f } ft = o and (1 _p) | at f _ i7 [o t } P ] f } _L f t = o ! (37) 

where 7 = —\J\- Obviously, a sufficient condition for a solution is 

af + [a,P]f = T Z\ and t f - i7[a t ,P]f = f Z 2 , (38) 

with some functions Zi(t,a,(v) and Z^it, a, cv). In the nonabelian case we choose Z\ = a and 
Z2 = while in the abelian case we take Z\ = z\ and Z2 = ^2 with z\i G C. In both these cases 
the compatibility conditions for the linear equations (^) lead to the equation 

[a\[a,P]] + [[a,P],[a\P]} = , (39) 

which is the field equation of the two-dimensional Euclidean sigma model [l7| . Since our P satisfies 
this equation, the equations (38) are compatible. After constructing a projector P by solving 
(38), one obtains a solution ip of the linear system ( |T3| ) with A and B obeying the field equations 
(15) and (|l~6|). Then one can choose ip, A, B as a new seed configuration and repeat the dressing 
transformation ip — > tp', again obtaining linear equations on some matrix T'. This iterative dressing 
procedure enables one to construct various solutions of the field equations. 
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4 Nonabelian multi-soliton configurations 



Equations. In order to generate some examples of nonabelian multi-solitons we specialize to U(2) 
and take as a seed configuration the simplest nontrivial solution of (p5|), viz. 



P 



1 



i 



TtT 



l+zz 
6 l+zz 



l+zz' 



1 



with 



T 



l+zz' 



and 



(40) 



We build the dressing factor x with a matrix f = ("), where u(t, z, z) and u(i, z, z) are functions 
to be determined. Substituting this ansatz for T into the first of eqs. (plf) with Zi = o, we obtain 



z, u 



z,v\ 



1 



LI + zz ' 
- 1 

2 L 1 + zz ' 



(u + zv) 



20 



l + zz 

. , 29 

(u + zi>) — z 



l + zz 



From these equations it follows that 

\z , zu — v] = 



zu - f(t,z) 



(41) 



(42) 



where f(t,z) is an arbitrary function of t and z (i.e. it does not depend on z). In a similar vein, 
the second of eqs. (|38|) with Z2 = reduces to 



1 



2, 



l + zz 



zv) 



Explicit solutions. It is not difficult to see that (f|l|) and ( fl3l) are solved by 
1 . 1 



1 + 



-zf 



-J 



with 



/ 



2i (t + 



l + zz"" 7 l + zz + 29' 

where h is an arbitrary meromorphic function of z (i.e. independent on t). Hence, we have 

2i 



1/ l + zz + 26 



(t + h(z)) 



(43) 



(44) 



(45) 



For locations given by t + h(z) = 0, we see that T = T and thus $ = (1 — 2P) 2 = 1 degenerates. 



One may obtain more general configurations by replacing T 



with T 



containing 



a meromorphic function g{z) and then solving (B8h again. In commutative limit this family of 
solutions coincides with the ones studied in |ji6|, [17] , |l8| 1 . For this reason we will not discuss them 
in detail. 

Depending on the explicit form of h(z), the solutions fliB] ) describe different kinds of multi- 
soliton configurations. Let us take, for example, h(z) = z q with q £ Z. Then, q < 1 leads to a time- 
dependent ring-like structure (for the energy density) while q > 2 creates a configuration consisting 
of q lumps. The latter simultaneously accelerate towards the origin (z = 0) of the noncommutative 
plane, scatter at an angle of n/q, and decelerate as they separate again (cf. |jl6|, [l?], |l8|j). 

Examples. To visualize the field configurations we employ the inverse Moyal-Weyl transforma- 
tion (10) and obtain 



and 



2i 

1 + zz + 



*(t + h(z)) (46) 
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with commutative coordinates t,z,z. Using these expressions we can calculate the projectors P± and 
P* as well as the field = (1— 2P*)*(1 — 2P*) and the energy density £+ for this family of solutions. 
Moreover, for large r 2 = zz the field and the energy density £+ will approach their commutative 



limits $ and £, respectively. Therefore, the asymptotic analysis of the papers [16, 17, 18] can be 
applied without alteration. 

The simplest scattering case occurs for h = z 2 . One can show that for large r the energy density 

peaks near the degeneration locus z 2 + i = 0. Hence, for a fixed large (positive or negative) time t, 
two lumps are centered at z = iy— t. Varying t we see the two lumps accelerating symmetrically 
towards each other along the x-axis, interacting at the origin near t=0 and decelerating to infinity 
along the y-axis. Thus, a head-on collision of these solitons results in 90° scattering. 

To exhibit a configuration with ring-like structure, one can take e.g. h = 0. In this case it turns 
out that the energy density 

is invariant under rotations around the origin. For fixed large time t, its degenerate maximum is 
found at r 2 ~ 2\t\. Varying t, this ring shrinks until t = and then expands again. One can also 
show that ± oo) — ► —1, i.e. no scattering occurs. 

In both these examples, noncommutativity is felt only inside a disc of radius R ~ \fd. For 
\t\ > then, the solitons agree qualitatively with their commutative cousins. 

5 Abelian multi-soliton configurations 

Explicit solution. The construction of abelian solitons is not achieved by simply taking the U(n) 
formalism to n = 1. The abelian case deserves a special treatment. For a seed solution, we take 
the simplest abelian static one-soliton configuration, made from the ket T=|0) and described by 
the projector 

P = |0)<0| (49) 



satisfying (25). The dressing factor \ is built from a ket T which, via (38), is subject to 



af— |0)(1| f = z\T and d t f - vy |1)(0| f = z 2 f . (50) 

We put z\ = Z2 = for simplicity^ Substituting a generic T = |0) + J2 n >i ^™(^)l n ) in t° these 
equations (the scale drops out), we learn that T n >2 = so that one is left with 

f = |0) + 2\(t)|l) =^ f = |0) + i 7 t|l) ■ (51) 

Therefore, we obtain 

P = T^fF = T -l I - I (|0)(0|+i 7 t|l)(0|-i 7 t|0)(l|+ 7 2 t 2 |l)(l|) . (52) 



1 More generally, Z\ and 22 characterize the locations of the solitons. 
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From the projectors ( |49|) and (|52|) we construct the field $ a la (p0|), 

* = 1 - rT ^2(7 2 i 2 |0)(0|+i 7 t|0)(l|+i 7 t|l)(0|+ 7 2 t 2 |l)(l|) . (53) 
It is easy to see that 

lim $ = 1 - 2(|0)<0| + |1)<1|) , (54) 
i.e. both limits coincide with the static two-soliton configuration. Therefore, no scattering occurs. 



Energy. The energy density for the solution (|53j ) is computed from (g), 
£ = %dt&dt$ + d g &dg$ + dg&d g & 

= \d t &d t § + ±[a\ &\[a\ $t]t + $t][ aj $t]t 

= 7 2 { ( IT ^ ¥ (|0)(0| + IDdl) + (i^lOXOI + ^(HXH + I2)(2|) 

- ( |2)(Q| + 10)(21 ) + (if 7 S) 2 (^ |1)(Q| ~ ^ |0)(11 + 12)(11 ~ |1)(2 } • (55) 

Consequently, the energy of this time-dependent abelian configuration is 

, f 4 4 7 2 t 2 4 7 2 t 2 

E = 2.8Tr£ = 2,8 7 j + + 

= (i^ 7 ^)l{ 1 + 2 ^ 2 + ^ 2 ) 2 } = 8 ^ = ' ( 56 ) 

where we used the definition 7 2 = 2/9. The value of 16n is twice the energy of a single static 
(rank-one) abelian soliton, in perfect agreement with the interpretation of (54). 

Inverse Moyal-Weyl map. The translation from operator- valued solutions &(t) to functions 
&±(t,z,z) solving the noncommutative field equations (|5|) is easily accomplished by the inverse 
Moyal-Weyl map (|i~0|). For the configuration (^) we obtain 

8 ' {2ft-i{z + z)}e-*fi , (57) 



1 



9 + 2t 2 

which produces the energy density 

16e-' r2 / e f „ „ 2r 2 /r 2 w 2t 2 \ 2 2/ 2 ,\ , 4 „ ,1 

By determining maximum values of this function one can see the ring-like structure of the solution 
and find the dependence of their size and shape on t. Note that these rings need not be symmetric. 
We see that our field ( |53| ) describes a two-soliton configuration, with both solitons sitting at the 
point x=y=0 and changing their shape with time. 

In order to construct more general time-dependent multi-solitons with a ring structure, one 
may start from a seed solution parametrized by (higher-rank) projectors Pk = |0)(0| + . . . + |A;)(A;| 
or P' k = Yli j=i \ a i) (ai\a ) ( a j\> wriere \cti) = e aiCtt_Qia |0) denote coherent states. In this context, 
we did not find abelian soliton scattering in the modified sigma model in 2+1 dimensions. To 
accomplish this, one should perhaps begin with time-dependent seed configurations or use a more 
general ansatz featuring ip with third-order poles in 



8 



6 Concluding remarks 



The present paper has applied a solution-generating technique (the dressing method) to the 
noncommutative modified U(n) sigma model in 2+1 dimensions. In this way we have constructed 
time-dependent classical field configurations which describe interacting noncommutative solitons. 
We have obtained ring-shaped bound states as well as nontrivial scattering solutions, albeit the 
latter ones only for the nonabelian case. It would be interesting to construct the Seiberg-Witten 
map for our soliton configurations as it was done e.g. in |l^, ^] for scalar solitons. Also, the brane 
interpretation of these solitons has not been completely clarified. 

Our construction is based on so-called dressing transformations, which map (old) solutions into 
(new) solutions. In fact, the set of all such transformations forms an infinite-dimensional group 



which acts on the space of all solutions of an integrable model (see e.g. [11, |T^] for discussion and 
references). Multi-solitons parametrize finite-dimensional orbits of this dressing group. 

As was explained in two previous papers || g] 3 the massless modes of the open N=2 string 
in a space-time filling brane with a constant NS S-field are described by noncommutative self-dual 
Yang-Mills (SDYM) theory in 2+2 dimensions. Upon reduction on a D2-brane worldvolume, the 
modified sigma model under consideration emerges. Therefore, the algebra of all (nonlocal) dressing 
symmetries of our field equations can be obtained (after noncommutative generalization) from the 
dressing symmetries of the SDYM equations (see e.g. |2^, ^4] and references therein) or from their 
stringy generalization fll^, Moreover, not only the modified sigma-model equations but a lot 
of other integrable equations in three and fewer dimensions derive from the SDYM equations by 



suitable reductions (see e.g. p6|, 27, 25, 29, pM and references therein). The investigation of the 



corresponding reductions of the noncommutative SDYM equations to noncommutative versions of 



KdV, nonlinear Schrodinger and some other integrable equations has just begun [31, [3^, 33, 34]. 
Their solitonic solutions are worth exploring. 
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